Large-scale tree-based unfitted finite elements for metal additive manufacturing

This thesis addresses large-scale numerical simulations of partial differential equations posed on evolving geometries. Our target application is the simulation of metal additive manufacturing (or 3D printing) with powder-bed fusion methods, such as Selective Laser Melting (SLM), Direct Metal Laser Sintering (DMLS) or Electron-Beam Melting (EBM). The simulation of metal additive manufacturing processes is a remarkable computational challenge, because processes are characterised by multiple scales in space and time and multiple complex physics that occur in intricate three-dimensional growing-in-time geometries. Only the synergy of advanced numerical algorithms and high-performance scientific computing tools can fully resolve, in the short run, the simulation needs in the area. The main goal of this Thesis is to design a a novel highly-scalable numerical framework with multi-resolution capability in arbitrarily complex evolving geometries. To this end, the framework is built by combining three computational tools: (1) parallel mesh generation and adaptation with forest-of-trees meshes, (2) robust unfitted finite element methods and (3) parallel finite element modelling of the geometry evolution in time. Our numerical research is driven by several limitations and open questions in the state-of-the-art of the three aforementioned areas, which are vital to achieve our main objective. All our developments are deployed with high-end distributed-memory implementations in the large-scale open-source software project FEMPAR. In considering our target application, (4) temporal and spatial model reduction strategies for thermal finite element models are investigated. They are coupled to our new large-scale computational framework to simplify optimisation of the manufacturing process. The contributions of this Thesis span the four ingredients above. Current understanding of (1) is substantially improved with rigorous proofs of the computational benefits of the 2:1 k-balance (ease of parallel implementation and high-scalability) and the minimum requirements a parallel tree-based mesh must fulfil to yield correct parallel finite element solvers atop them. Concerning (2), a robust, optimal and scalable formulation of the aggregated unfitted finite element method is proposed on parallel tree-based meshes for elliptic problems with unfitted external contour or unfitted interfaces. To the author’s best knowledge, this marks the first time techniques (1) and (2) are brought together. After enhancing (1)+(2) with a novel parallel approach for (3), the resulting framework is able to mitigate a major performance bottleneck in large-scale simulations of metal additive manufacturing processes by powder-bed fusion: scalable adaptive (re)meshing in arbitrarily complex geometries that grow in time. Along the development of this Thesis, our application problem (4) is investigated in two joint collaborations with the Monash Centre for Additive Manufacturing and Monash University in Melbourne, Australia. The first contribution is an experimentally-supported thorough numerical assessment of time-lumping methods, the second one is a novel experimentally-validated formulation of a new physics-based thermal contact model, accounting for thermal inertia and suitable for model localisation, the so-called virtual domain approximation. By efficiently exploiting high-performance computing resources, our new computational framework enables large-scale finite element analysis of metal additive manufacturing processes, with increased fidelity of predictions and dramatical reductions of computing times. It can also be combined with the proposed model reductions for fast thermal optimisation of the manufacturing process. These tools open the path to accelerate the understanding of the process-to-performance link and digital product design and certification in metal additive manufacturing, two milestones that are vital to exploit the technology for mass-production.Aquesta tesi tracta la simulació a gran escala d'equacions en derivades parcials sobre geometries variables. L'aplicació principal és la simulació de procesos de fabricació additiva (o impressió 3D) amb metalls i per mètodes de fusió de llit de pols, com ara Selective Laser Melting (SLM), Direct Metal Laser Sintering (DMLS) o Electron-Beam Melting (EBM). La simulació d'aquests processos és un repte computacional excepcional, perquè els processos estan caracteritzats per múltiples escales espaitemporals i múltiples físiques que tenen lloc sobre geometries tridimensionals complicades que creixen en el temps. La sinèrgia entre algorismes numèrics avançats i eines de computació científica d'alt rendiment és la única via per resoldre completament i a curt termini les necessitats en simulació d'aquesta àrea. El principal objectiu d'aquesta tesi és dissenyar un nou marc numèric escalable de simulació amb capacitat de multiresolució en geometries complexes i variables. El nou marc es construeix unint tres eines computacionals: (1) mallat paral·lel i adaptatiu amb malles de boscs d'arbre, (2) mètodes d'elements finits immersos robustos i (3) modelització en paral·lel amb elements finits de geometries que creixen en el temps. Algunes limitacions i problemes oberts en l'estat de l'art, que són claus per aconseguir el nostre objectiu, guien la nostra recerca. Tots els desenvolupaments s'implementen en arquitectures de memòria distribuïda amb el programari d'accés obert FEMPAR. Quant al problema d'aplicació, (4) s'investiguen models reduïts en espai i temps per models tèrmics del procés. Aquests models reduïts s'acoplen al nostre marc computacional per simplificar l'optimització del procés. Les contribucions d'aquesta tesi abasten els quatre punts de dalt. L'estat de l'art de (1) es millora substancialment amb proves riguroses dels beneficis computacionals del 2:1 balancejat (fàcil paral·lelització i alta escalabilitat), així com dels requisits mínims que aquest tipus de mallat han de complir per garantir que els espais d'elements finits que s'hi defineixin estiguin ben posats. Quant a (2), s'ha formulat un mètode robust, òptim i escalable per agregació per problemes el·líptics amb contorn o interface immerses. Després d'augmentar (1)+(2) amb un nova estratègia paral·lela per (3), el marc de simulació resultant mitiga de manera efectiva el principal coll d'ampolla en la simulació de processos de fabricació additiva en llits de pols de metall: adaptivitat i remallat escalable en geometries complexes que creixen en el temps. Durant el desenvolupament de la tesi, es col·labora amb el Monash Centre for Additive Manufacturing i la Universitat de Monash de Melbourne, Austràlia, per investigar el problema d'aplicació. En primer lloc, es fa una anàlisi experimental i numèrica exhaustiva dels mètodes d'aggregació temporal. En segon lloc, es proposa i valida experimental una nova formulació de contacte tèrmic que té en compte la inèrcia tèrmica i és adequat per a localitzar el model, l'anomenada aproximació per dominis virtuals. Mitjançant l'ús eficient de recursos computacionals d'alt rendiment, el nostre nou marc computacional fa possible l'anàlisi d'elements finits a gran escala dels processos de fabricació additiva amb metalls, amb augment de la fidelitat de les prediccions i reduccions significatives de temps de computació. Així mateix, es pot combinar amb els models reduïts que es proposen per l'optimització tèrmica del procés de fabricació. Aquestes eines contribueixen a accelerar la comprensió del lligam procés-rendiment i la digitalització del disseny i certificació de productes en fabricació additiva per metalls, dues fites crucials per explotar la tecnologia en producció en massa.Postprint (published version

Treating inextensibility constraints in hyperelastic materials with inequality level sets

Elastomers are viscoelastic polymers with low Young's modulus and high failure strain that are used in many c ivil e ngineering applications, including bridge bearings, seismic isolators for buil dings and resilient rail wheels. Their constitutive behaviour i s characterized by a nonlinear stress - strain relation with an extensibility limit . This contrasts with materials that have instead a limit on the tensile stresses, such as mild steel. This MSc thesis is concerned with the numerical modeling of elastomers. T his involves dealing with a medium with two phases: a constrained region, wher e the particles have reached their maxi mum allowable deformation, and a free region, where the inextensibility constraint is still inactive. Moreover, one can think of an interf ace splitting the two phases of the medium. If the focus is put in obtaining methods to locate and evolve such interface, then a two phase medium with a moving interface is considered. From the mathematical point of view, this is a constrained minimization problem. One of the strategies to solve it is to turn the minimization problem into a shape equilibrium one. This approach has b een successfully employed for an interface location problem in small strains and serves as the starting point of this work. T hu s, t he main purpose of this thesis is to extend this formulation to a large strains interface locating and evolving scenario. A first analysis of the problem reveals two sources of nonlinearity : the inextensibility constraint and the kinematics in large st rains. A simple but thorough one - dimensional study of the problem is then developed to find methods to sort out both nonlinearities. Following this, explicit iterative schemes to locate and evolve one or multiple interfaces are straightforwardly obtained i n 1D linear elasticity. However, the same ideas applied to a simple St.Venant - Kirchhoff hyperelasticity model , evidences that even very simple 1D problems become rather complex and cannot be solved as directly and explicit as before. Numerical examples are provided throughout this analysis and they are also useful to conclude that both locating and evolving the interface can be essentially seen as the same problem, but with different driving effects. After that, an extension of the one - dimensional schemes t o two or more dimensions is explored . Although the same ideas can be applied, more sophisticated modeling tools are required, namely, the X - FEM and Level set methods, the shape sensitivity analysis and the Arbitrary Lagrangian - Eulerian methods . A complemen tary numerical implementation of the pr oposed strategy is to show its computational benefits. In particular, a combination of the three previous techniques shall make unnecessary a stepwise update of the Level set. T he work presented here may not be limite d to this particular case and be r elevant to other engineering problems involving moving interfaces and boundaries , such as plasticity analysis or the saturation of a porous medium

A generic finite element framework on parallel tree-based adaptive meshes

We present highly scalable parallel distributed-memory algorithms and associated data structures for a generic finite element framework that supports h-adaptivity on computational domains represented as multiple connected adaptive trees—forest-of-trees—, thus providing multi-scale resolution on problems governed by partial differential equations.The framework is grounded on a rich representation of the adaptive mesh suitable for generic finite elements that is built on top of a low-level, light-weight forest-oftrees data structure handled by a specialized, highly parallel adaptive meshing engine. Along the way, we have identified the requirements that the forest-of-trees layer must fulfill to be coupled into our framework. Essentially, it must be able to describe neighboring relationships between cells in the adapted mesh (apart from hierarchical relationships) across the lower-dimensional objects at the boundary of the cells. Atop this two-layered mesh representation, we build the rest of data structures required for the numerical integration and assembly of the discrete system of linear equations.We consider algorithms that are suitable for both subassembled and fully-assembled distributed data layouts of linear system matrices. The proposed framework has been implemented within the FEMPAR scientific software library, using p4est as a practical forest-of-octrees demonstrator. A comprehensive strong scaling study of this implementation when applied to Poisson and Maxwell problems reveals remarkable scalability up to 32.2K CPU cores and 482.2M degrees of freedom. Besides, the implementation in FEMPAR of the proposed approach is up to 2.6 and 3.4 times faster than the state-of-the-art deal.II finite element software in the h-adaptive approximation of a Poisson problem with firstand second-order Lagrangian finite elements, respectively (excluding the linear solver step from the comparison)

Numerical modelling of heat transfer and experimental validation in Powder-Bed Fusion with the Virtual Domain Approximation

Among metal additive manufacturing technologies, powder-bed fusion features very thin layers and rapid solidification rates, leading to long build jobs and a highly localized process. Many efforts are being devoted to accelerate simulation times for practical industrial applications. The new approach suggested here, the virtual domain approximation, is a physics-based rationale for spatial reduction of the domain in the thermal finite-element analysis at the part scale. Computational experiments address, among others, validation against a large physical experiment of 17.5 $\mathrm{[cm^3]}$ of deposited volume in 647 layers. For fast and automatic parameter estimation at such level of complexity, a high-performance computing framework is employed. It couples FEMPAR-AM, a specialized parallel finite-element software, with Dakota, for the parametric exploration. Compared to previous state-of-the-art, this formulation provides higher accuracy at the same computational cost. This sets the path to a fully virtualized model, considering an upwards-moving domain covering the last printed layers

Distúrbios metabólicos por manejo alimentar inadequado em ruminantes: novos conceitos

In the “top” cattle farming, the metabolic problems are frequent due to the necessity of getting an animal that presents a high performance. One of the main problems found in this type of activity is the extreme concentration of food supply. Among the main metabolic disorders are ruminal acidosis, bloat, ketosis and hepatic lipidosis. The metabolic disorders’ damages are rarely visible because, most of times, these illnesses don’t show clinical signals. There are ways to manipulate the rumen to minimize the metabolic disorders that can be at the animal level, the diet or rumen microorganisms, such as the use of residues and by-products from the agro-industry in the animal feeding. More researches are necessary to completely understand the influence of the nutrition handling on the metabolic disorders.En las llamadas ganaderías “elite” los problemas metabólicos son frecuentes debido a la necesidad de obtener un elevado desempeño animal. Uno de los principales problemas observados en esos sistemas es la elevada oferta de alimentos concentrados. Dentro de los principales desórdenes metabólicos se destacan la acidosis ruminal, timpanismo, cetosis e hígado graso. Los perjuicios ocasionados por estos problemas no son visibles, ya que muchas veces no aparecen síntomas clínicos. Existen formas de manipular el rumen para evitar estos trastornos, a través de la dieta, del animal o de los microorganismos ruminales, como puede ser el uso de ciertos alimentos. Son necesarias más investigaciones sobre la influencia del manejo alimenticio para evitar la incidencia de estos disturbios y poder entenderlos mejor

AS CONTRIBUIÇÕES DA QUARTA REVOLUÇÃO INDUSTRIAL PARA ECOINOVAÇÕES: UMA REVISÃO BIBLIOMÉTRICA

A partir desenvolvimento dos setores da tecnológicos, avista-se inovações que têm se mostrado satisfatórias, sobretudo diante dos efeitos da urbanização. No entanto, cada vez mais se vê a necessidade da modernização tecnológica focada em sustentabilidade, de modo a alcançar equilíbrio entre o ambiental, o social e o econômico. Diante disso, levanta-se a questão de se saber qual o verdadeiro impacto das preocupações ambientais na inovação desenvolvida e adotada pelos diversos agentes econômicos. O objetivo deste artigo, portanto, é identificar como as tecnologias trazidas pela Quarta Revolução Industrial podem contribuir para a adoção de estratégias de ecoinovação.[1] Center for Sustainable Development, University of Southern Santa Catarina (UNISUL), e-mail:

Characterization of Spontaneous Bone Marrow Recovery after Sublethal Total Body Irradiation: Importance of the Osteoblastic/Adipocytic Balance

Many studies have already examined the hematopoietic recovery after irradiation but paid with very little attention to the bone marrow microenvironment. Nonetheless previous studies in a murine model of reversible radio-induced bone marrow aplasia have shown a significant increase in alkaline phosphatase activity (ALP) prior to hematopoietic regeneration. This increase in ALP activity was not due to cell proliferation but could be attributed to modifications of the properties of mesenchymal stem cells (MSC). We thus undertook a study to assess the kinetics of the evolution of MSC correlated to their hematopoietic supportive capacities in mice treated with sub lethal total body irradiation. In our study, colony-forming units – fibroblasts (CFU-Fs) assay showed a significant MSC rate increase in irradiated bone marrows. CFU-Fs colonies still possessed differentiation capacities of MSC but colonies from mice sacrificed 3 days after irradiation displayed high rates of ALP activity and a transient increase in osteoblastic markers expression while pparγ and neuropilin-1 decreased. Hematopoietic supportive capacities of CFU-Fs were also modified: as compared to controls, irradiated CFU-Fs significantly increased the proliferation rate of hematopoietic precursors and accelerated the differentiation toward the granulocytic lineage. Our data provide the first evidence of the key role exerted by the balance between osteoblasts and adipocytes in spontaneous bone marrow regeneration. First, (pre)osteoblast differentiation from MSC stimulated hematopoietic precursor's proliferation and granulopoietic regeneration. Then, in a second time (pre)osteoblasts progressively disappeared in favour of adipocytic cells which down regulated the proliferation and granulocytic differentiation and then contributed to a return to pre-irradiation conditions